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Introduction. An analytical grasp of the thermody- 
namics of SU(3) gauge theory is needed to interprete 
ultrarelativistic heavy ion collisions. Thermal pertur- 
bation theory (TPT) by itself is unable to provide a 
framework for this. Nonperturbative lattice input is 
needed for a prediction of the thermal pressure . Lat- 
tice experiments are a nonperturbative, first-principle ap- 
proach. Although numerical in nature, they have boosted 
our knowledge about relevant degrees of freedom and 
their condensation. It is by now an established result 
that condensed center vortices are the relevant degrees 
of freedom for confinement at temperature T below Tc 
0- The purpose of this Letter is to develop a nonper- 
turbative approach to thermal SU(iV) thermodynamics 
based on thermal quasiparticle excitations (QPEs) and 
the condensation of the solitonic objects caloron, mag- 
netic monopole and center vortex (pointlike in 4D, 3D, 
and 2D Q). The latter processes successively reduces 
the asymptotic gauge symmetry SU(N) according to the 
pattern SU(N)-^U(1)~-^^Zn ^nothing. The first two 
symmetry reductions are accompanied by a 2°'' (possi- 
bly occuring close to the Planck mass) and a very mild 
l""* order thermal phase transition, respectively. Order 
parameters for these transitions are the masses of gauge 
bosons associated with coset generators and the Cartan 
subalgebra, respectively. The final transition towards 
spontaneously broken, local symmetry is 1°' order. 
The order parameter is the energy density p. A single 
mass scale A^ or Am or Ac, determined by a boundary 
condition, introduces a gauge invariant scale separation 
in all phases. 

Analysis. At temperatures very much higher than the 
dynamic scale Aym(N) TPT, when applied to the calcu- 
lation of thermodynamic potentials, shows convergence 
in its perturbative orders. 

Electric phase. At a critical temperature T^, much 
larger than Aym(N), a thermal 2'"* order phase transi- 
tion occurs at which the system condenses a compos- 
ite, adjoint Higgs field (j> reducing the gauge symme- 
try maximally, SU(N)— > U(l)'^ (lowering T across the 
transition, the order parameter gauge boson mass con- 
tinuously increases starting from a very small screening 
mass ). Relevant degrees of freedom driving this con- 
densation are calorons. At high temperature an isolated 
caloron is a pointlike object with a nontrivial holonomy 
and topological charge one in 4D Euclidean spacetime. 
A caloron possesses N Bogomoln'yi-Prasad-Sommerfield 
(BPS) magnetic monopole constituents 0, 0| related by 



center transformations. Its action is S'^^i = 87r^/e^ j^, 
where e denotes the gauge coupling. The caloron is a 
BPS saturated solution to the SU(N) Yang-Mills equa- 
tion of motion , its 3D energy density vanishes. It is a co- 
herent thermal state in 4D. The picture that caloron con- 
densation causes the growth of gauge boson masses with 
decreasing T is self-consistent if S^^i rapidly decreases 
below by a large rise of the gauge coupling e. This 
indeed happens, see Fig. la. An indirect argument in fa- 
vor of caloron condensation relies on the lattice result of a 
vanishing 4D topological charge density x for T > im- 
plying that calorons do not populate the ground state as 
small-size objects. The caloron radius and e are strictly 
monotonously related. We conclude that above T^. e is 
large and the caloron action small. After caloron con- 
densation has taken place, the proposed effective action 
is 

Se^J^ dT j d\ Q tr G,,,G^, -I- trV^c^V^c^ + Ve(0))(1) 
where Ve denotes the potential for the field 



variant derivative is defined as 'D^cj) — 



The co- 

and the field strength as G^i/ — G°^<°, where G°j, = 
d^Al - d^A^ ~ ef^^A^^Al and tv^tH'' = l/2(5''^ In a 
suitable gauge the field describes the coherent ground 
state of the thermal system. As a consequence, its gauge 
invariant modulus is independent of spatial coordi- 
nates. A dependence oi\(j)\ on r may only enter through 
an adiabatically slow r dependence of the temperature T . 
In the absence of gauge fields in the 3D energy density 
of (j) is exactly zero since it is composed of BPS saturated 
calorons. This condition is satisfied for the configuration 
(j) if and only if it is BPS saturated itself. In other words: 
we search a potential Ve such that (i) a BPS saturated 
and spatially homogeneous solution to the classical equa- 
tion of motion, derived from Q in the absence of gauge 
fields, exist, (ii) the modulus |0| of this configuration is 
independent of r, and (iii) quantum and thermal fluctu- 
ations of (j) are negligible. We proceed by discussing the 
case of even N. The other case is more complicated, an 
analysis is given in Here results for N=3 are quoted 
in square brackets if appropriate. We work in a gauge 
where 4> = diag(0i, 02, • • • , ^N/2) is SU(2) block diagonal 
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zA|;diag(Ai0i/|(/.i|^-- - ,Ai0«/2/|0N/2r). (2) 



The SU(2) modulus is defined as = 1/2 tr2 0; (/ = 
1, • • • ,N/2). The gauge invariant potential Ve is given 



as Vsicl)) = Altr^ 



^. The periodic solution to the 
= ve (Ve = tif^v'^) with minimal 



BPS equation = ve {Ve 
potential and maximal SU(iV) symmetry breaking is in 
SU(2) decomposition given as 



Mr) 



K%/{2ttTI) A3 exp(-27rfrZAir) . 



(3) 



hermitian matrix. In uni- 
^n), there is no r de- 



Obviously, </) is a traceless, 
tary gauge, cj) = diag(0i , </)2 , 

pendence (no r dependence occurs in trjv</)^). From ^ 
and ISI) we obtain the following ratios: d^- .Ve/T'^ = 

127r2p and |Ve/|0/P = 3/^A|; where the dimension- 
less temperature is defined as A_b = 2'kT / Ke ^ 1- 
As a consequence, thermal fluctuations of cf) are negli- 
gibly small and quantum fluctuations do not exist (the 
mass of (pi excitations is larger than the composite- 
ness scale \4>i\)- A pure-gauge solution to the equation 
of motion 'Df_iGf^i, = 2ie[(f),'Di,(j3\ in the background of 
® is A^a = (7r/e)T(5^odiag(Ai,2Ai,-- - ,N/2Ai). We 
have Vfj^cj) — for A'lf entering in P^. As a conse- 
quence, the (j) kinetic term in the action ^ vanishes on 
(j) and A'lf , the ground-state energy density is lifted to 
Vj;((/,) = 7r/2A|rN(N-|-2) [T4;(0) = 47rA|r] by caloron 
interaction. The field (j) back-reacts by emitting and ab- 
sorbing the thermal quasiparticle excitations (QPEs) de- 
scribed by fluctuations a^. This makes sense since (j) is 
composed of coherent thermal states - calorons. Substi- 



tuting A^ 



into the action its form is 



retained as an action for the fluctuations (the ki- 
netic term induces mass terms for a^^ only in unitary 
gauge, where J^^ = 0. QPEs are related to fluctua- 
tions only in this gauge. Unitary gauge is reached 



by a gauge transformation Ay^ Of^^fi -I- i/e((9^ri^)ri 
(fi = e*^) involving a nonperiodic 9 with SU(2) decompo- 
sition 61 = —ttXiTIt. Because the matter field cj) does not 
fluctuate the physics is inert under such a gauge trans- 
formation: Hosotani's mechanism |3| does not take place. 
Out of the independent modes N-1 modes reside in 
the Cartan subalgebra and remain massless on tree level 
(TLM modes). In unitary gauge the N(N-l) independent 
off-diagonal modes, defined w.r.t. the SU(N) generators 
tii = 1/2 iS'rS'I+S'X) andt^f = (SiSi-SiS^) {I = 
1, • • • ,N) {J > I), are heavy on tree level (TLH modes). 
Their mass spectrum is obtained from the (j) kinetic term 
in (HJ as = -26^ tr t"] (^ = {t",t"}). We 



have ; 



r)2 for both t" and t" . Since TLH 



and TLM modes are QPEs radiative corrections are iocaJ 
in unitary gauge, unitarity is trivially satisfied. A radia- 
tive correction Amj j to j and a radiatively induced 
mass-squared for a TLM mode rri^j^j^j are determined by 
a single 4-gauge boson tadpole with a TLH or a TLM 
mode and a TLH mode running in the loop, respectively. 
Since a loop integral is cut off at the compositeness scale 
\(f>\, the expansion of a radiative correction is obtained as 

'"e~'^-\ ) and 

•). ForN=2 we 



Am2^ = |0|V(8^2)(c-\«e2 + c, 



„TLH„0 
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I0l7(8^^)(cj^ 



have Cq^" = 1/8. Radiative corrections AVg to the po- 
tential Ve are induced by one-loop and two-loop bubbles. 




FIG. 1: The evolution of the gauge coupling in the electric (a) 
and magnetic phase (b) for N=2 (thick grey lines), N=3 (thick 
black lines), N=4 (dashed lines), and N=10 (thin solid lines). 
The initial conditions in (b) are obtained by matching the 
pressure at the phase boundary. In both cases the gauge cou- 
plings diverge logarithmically, e,g ^ — log(AB,A/ — ^'e^m)^ b± 
Xe = {9.92,4.79,8.6,7.87} and A^^ = {3.55,3.36,4.48^4-64}, 
respectively. This is not captured by Mathematica numerics. 



In the latter case two one-loop bubbles are connected by 
a 4-gauge boson vertex. AVe is tiny. To lowest oder in 
l/e2 we have /^Ve/Ve < (N^ - l)/(16^2)(|^|/^^)6 ^ i. 
see below. We proceed by considering tree-level masses 
only. A useful quantity is — Cka (mass/T of a TLH 
mode) where a = e^J K%,/{2nT'^) (or a = g K\J {2ttT^) , 
see below). The value of the dimcnsionless constant 
Ck (k = I,-- - ,N(N-1)) derives from the TLH mass 
spectrum. TLH modes are thermal QPEs - they ac- 
quire mass and structure by interacting with the ground 
state of the system. Thus the thermodynamic relation 
Pe = — Pe between the total energy density pe{T) 

and pressure Pe{T) does not follow from the partition 
function associated with the action ^ for the fiuctua- 
tions flp . This problem is resolved by imposing the condi- 
tion of minimal thermodynamic self-consistency daP ~ 
jSj. An evolution equation follows as 



da^E — — 



24 Al 



(27r)6N(N+2; 



N(N-l) 

■E 



clD{ak) , 



(4) 



where D (a) = dxx^ l{\Jx^ -I- a^(exp(-\/a;^ -I- a^) — 1)). 
The right-hand side of Q is negative definite: the 
function A£;(a) can be inverted to a(\E\ As a con- 
sequence, the gauge coupling is obtained as e(A£;) = 

a(A£;)A'^''^/(27r). To compute e(A£;) numerically, subject 
to the initial condition a{\^ — 1000) = 0, we have de- 
vised a Mathematica program. The result is shown in 
Fig. la. We have checked numerically that the value of 
\% does not depend on the initial condition a(A|^) = 
within a wide range of A^ values. The plateau values of 
e are e ~ {17.15,8,14,12} for N=2,3,4,10, respectively. 
We conclude that the caloron action S^^i — 87r^/e^ is 
very small throughout the entire electric phase. The as- 
sumed condensation of calorons, driving the condensa- 
tion of 0, is self-consistent by a large rise of e across 
the phase boundary. There is a regime in A^; where e is 
large but much smaller than its value close to the critical 
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values A|; and A|^. Calorons are of finite size but still 
large enough not to generate a topological susceptibility 
on state-of-the-art lattices. Constituent BPS magnetic 
monopoles become visible [Bj. Calorons start to scat- 
ter elastically (SU(2) calorons of the same embedding 
in SU(N) go through each other), scatter inelastically 
(calorons of different embedding approach), and annihi- 
late (caloron meets its anticaloron) . No monopoles, sin- 
gle monopoles, and monopole-antimonopole pairs arise, 
respectively. Monopoles are slowed down by surrounding 
calorons. A pair separates into stable particles connected 
by a magnetic flux line, a long-lived dipole forms. If a sin- 
gle monopole is produced then it is unstable (magnetic 
charge is conserved!) unless it connects with its anti- 
monopole, produced in an independent collision. There 
are N species of monopoles equal to the maximal num- 
ber of monopole constituents a free-of-magnetic-charge 
SU(N) caloron can have. How do we understand this 
in the effective electric theory? In the gauge, where (p 
is SU(2) block diagonal, the temporal winding of the 
l*^ block is complemented by spatial winding at isolated 
points in 3D space. By a large, r dependent gauge trans- 
formation the monopole's asymptotic SU(2) Higgs field is 
rotated to spatial constancy and into the direction given 
by the temporal winding of the unperturbed block (j)i. 
Its Dirac string rotates in Euclidean time. Since this 
monopole is stable, a correlated antimonopolc must exist. 
We arrive at a dipole rotating about its center of mass at 
an angular frequency 2nlT. This rotation is an artifact of 
our gauge choice. Rotating the dipole to unitary gauge by 
the above gauge function 0;, we arrive at a (quasi)static 
dipole. Recall that the fluctuations are thermal QPEs 
only in unitary gauge. There are isolated coincidence 
points (CPs) in time where the lower right (upper left) 
corner tpi (ipi+i) of the l'" {{I + 1)*") SU(2) block (now 
(Z = l,---,N/2 — 1)) together with the number zero of its 
right-hand (left-hand) neighbour are proportional to the 
generator A3 . Coincidence also happens between the first 
and last diagonal entry in (j). At a CP, spatial winding 
may take place at isolated points in 3D space. Monopoles 
of this type are unstable particles. Coincidence also hap- 
pens between nonadjacent SU(2) blocks, see Fig. 2. The 
associated monopoles are, however, not independent. We 
obtain N/2-|-N/2=N monopole species - in accord with 
the above counting [The field (f) winds with winding num- 
ber one in each of the independent SU(2) subalgebras 
for half the time. The matching happens at the CPs 
TCP — 0, l/(2r) where 4>{tcp ± 0) add up to an element 
of the third, dependent SU(2) algebra. At these CPs we 
may have an unstable monopole.]. Monopoles have res- 
onant excitations and are coherent thermal states in a 
thermal environment [Tol |. 

Magnetic phase. The mass Mi of a stable (unstable) 
BPS monopole is given as Mi ~ 8(l/-\/2)7r|0/ |/e. It van- 
ishes at A|;. As a consequence. Mi and the monopole 
action Mi/T vanish. TLH modes decouple at X'^. Mag- 
netic monopoles condense in a 2°'' order phase transition 
(radiative corrections in the electric phase lift this to a 
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FIG. 2: Sum of (p2 and ip^ (bold curve) and sum of 992 and 
ifs as a function of r. A zero in each curve is associated with 
an unstable monopole or its unstable antimonopole. 



very weak 1°' order transition). The order parameter 
is the TLM mass. It is continuous across the transi- 
tion. U'*"^ is spontaneously broken by the condensates 
of stable monopoles described by the complex fields t/f/c 
(fc — I,-- - ,N/2[2]). These fields are winding in time 
with integer winding number. The center symmetry Zjv 
is unbroken. It is represented by a local permutation 
symmetry (j)k ~^ (j)k+j (i(x) = 1, • • • , N - 1; for fc > N/2 
it maps zero onto zero) under which the ground state 
is inert. The potential for the complex field 0^ is de- 
fined as V^j^ = ' where w^^^ = iK\J(j3k- The 
ground state of the system possesses an energy den- 
sity 1/2 Fm = (^/8)TA3,N(N+2) [1/214/ = ^TAf,] 
(canonical definition of 0-kinetic term). It is constructed 
in close analogy to the electric phase. Again, it is as- 
sured that there are no thermal and quantum fluctu- 
ations of the fields 4>k- In unitary gauge, 4>k = j^fej, 
we define QPEs as thermal fluctuations of the Abelian 
gauge bosons. These QPEs are all massive away from 
the phase boundary. Their mass spectrum is ruk = g\4>k\^ 
where g denotes the magnetic gauge coupling. At the 
electric-magnetic phase boundary, we have g = 2tt /e = 
0. There are no radiative corrections to the gauge bo- 
son masses in the magnetic phase. Minimal thermo- 
dynamic self-consistency implies an evolution equation 
for 5: da\M = -96Al^a/((27r)6N(N-f2))^^f,4D(afe) 
[da\M = -l'2\\,iaD{a)/{2TTf]. The evolution of g with 
temperature is shown in Fig. lb. Monopole condensa- 
tion is driven by a large rise in the gauge coupling g 
across the phase boundary. The size of a BPS magnetic 
monopole is given by a constant C^^ (x g 11]. Accord- 
ing to Fig. lb the size of a magnetic monopole increases. 
Isolated dipoles and unstable single monopoles 'feel' each 
other close to the boundary to the center phase, where g 
diverges logarithmically. The former form stable strings 
of dipoles. The latter form stable dipoles first and then 
stable strings of dipoles. Strings of dipoles are interpreted 
as fat center vortices 13] . Monopole constituents of a vor- 
tex are coherent thermal states and so is a vortex. 
Center phase. Although a center vortex and an 
Abrikosov-Nielsen-Oleson (ANO) vortex are not the 
same solitonic object, the Euclidean action and the en- 
ergy of the former can be estimated by the latter. 
For a typical ANO vortex (length ~ width) they are 
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and E,„, cx l/g^ [H. Since g 



given as 

diverges logarithmically at center vortices condense 
and Abelian gauge bosons decouple. The equation of 
state (EOS) at X%i is p = —P. The condensed phase 
is described by complex and 3D local fields $„ , n = 
f , • • • N. These fields are, up to a dimensionful (and T 
dependent) normalization, the condensates of spherical, 
spatial Wilson loops of diameter comparable to the core- 
sizes of the respective center vortices. The (local!) 



n+p 



symmetry acts as 
(p = p(2^) = Ij ■ ■ ■ T^)- Thus it is spontaneously broken. 
The thermodynamics of the fields (n = 1, • • • , N) is 

determined by symmetric potentials Vf 



(")- 



where 



„(") 



ate minima of the potential VX at 



-(") (") 
The N degener- 

'i'nl = are ex- 
actiy at zero. At A^^ a jump of the order parameter 
p to p = takes place. The latent heat ec is given as 
ec - (N(N+2)/16)A$;^Al^ [ec = (1/2)A^mA4^]. The pres- 
sure P in the center phase starts at the negative value 
—2ec and relaxes rapidly to zero. If a matching to the 
magnetic phase would take place in thermal equilibrium 
no QPEs would exist below A^^ since the potential is well 
approximated by its pure pole term, and the ussual cal- 
culation for d^^ i^c"'' applies. This situation holds for 
|$„| sufficiently smaller than Ac or for very large N. The 
solutions to the BPS equations would then be identical 
to the ones in the magnetic phase. Deviations from ther- 
mal equilibrium at the center transition manifest them- 
selves in a r dependence of see Since is 
a local symmetry no domain walls are generated in the 
transition. Heavy 'glueballs' are produced by a rapid re- 
heating, see 20] for the interpretation of these 'hadrons'. 
Matching the phases. P and T are continuous across 
a thermal phase boundary. In the electric and mag- 
netic phase the ratio of gauge and kinetic terms de- 
fines the mass of the gauge-field spectrum at a given T. 
A re-definition of this ratio re-defines the scales A^ and 
Am, respectively. Using the canonical definition (also for 
the center phase) and disregarding a polarization mis- 
match for Abelian gauge bosons at the electric-magnetic 
phase boundary, this relates the scales Ae, Am , Ac as 
Ae = (1/4)1/3Am = (l/4)i/3Ac [Ae = (1/2)1/3Am = 
(1/3)1/3 Ac] on tree-level. 

Energy density p and pressure P. In Fig. 3 the T 
dependence of p and P is shown. The results for other 
thermodynamic potentials and for the critical exponents 
of the perturbative-electric and electric-magnetic transi- 
tion will be published elsewhere. The value of p at A^^ is 
a measure of the latent heat released in the 1°' order-like 
transition following the onset of center-vortex condensa- 
tion. 

Deviation from thermal equilibrium (DTE) at \^ : A 
measure for this deviation is the relative deviation of 
the full potential Vc from its pure pole term evaluated 
on the lowest winding mode taken at Am = X%i- For 
N = 2,3,4,10 we respectively obtain: 70%, 41%, 23%, 
and 1%. The large DTE for N=2 may explain the appar- 



l^as /T ^ = 5.48 
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FIG. 3: The normalized energy density p/T* and the normal- 
ized pressure P/T^ as functions of Ab. The horizontal lines 
indicate the respective asymptotic value (taken at Ab = 100). 
The discontinuities are due to the jump in gauge-boson po- 
larizations across the electric-magnetic transition, this effect 
disappears beyond tree-level (Abelian gauge bosons in the 
electric phase are then massive). Notice the brevity of the 
magnetic phase. 



ent 2"'' order deconfinement transition seen on the lattice. 
It may be an artifact of imposing thermal equilibrium in 
the lattice simulation. For N=3 the DTE is still large. It 
is, however, conspicuous that a high-statistics determi- 
nation of P, not relying on the method of (3 integration, 
obtains a negative pressure slightly above Tc 14]. The 
/3-integration method may be biased by a prescribed in- 
tegration constant. For N=4 and in particular for N=10 
the situation is much better for lattice simulations. It is 
extremly likely that future lattice simulations with N>3 
will unambiguously identify negative pressure close to A^ 
(our confidence regarding the truth of this prediction de- 
rives from the results in 20]). For N=2,3,4,10 we predict 
negative pressure of magnitude as given in Fig. 3. 
Comparison with the lattice. For comparison with 



iin: (Pi„t/p)j 

[13 Gil 113 



For 



to correspond to A' 
(T= 1.535; 4.878; 12.925 
(18^ X 5; 10^ X 3; 10^ x 5 lat 



lattice data we take 
N = 2 we obtain: 

0.31; 1.26; 1.29 ti^ttiaiii* - 

tices), {Piat/P){T = 4.878 Tc) = 1.26 O (10^ x 3 
lattice). For T = 12.195 Tc ^'lat is calculated using 
Ai„t = Put - 3Pi„t and piat- For iV = 3 we obtain: 
Piat/p{T = 5r,) = 0.88, PiatlP{T = 5T,) = 0.74 [3 
(32^ X 8 lattice). Deviations may be explained by the 
tree-level treatment of gauge boson masses in the elec- 
tric phase and a mismatch with the lattice scale [Tgl ]. 
Summary. We have developed a nonperturbative, in- 
ductive approach to thermal SU(N) Yang-Mills theory. 
The numerical results for pressure and energy density are 
promising. A more precise computation will aim at a pre- 
cise match of the lattice scale and the scale Ac including 
radiative corrections in the electric phase. This predicts 
the critical temperature Tp in the UV regime (possibly 
close to the Planck mass) . The detection of negative pres- 
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sure at T close to Tc has far-reaching consequences for 
our understanding of accelerated cosmological expansion 
and the creation of a lepton asymmetry ;20|. QCD ther- 
modynamics will be described by an extension of the ap- 
proach which includes dynamical, fundamental fermions, 
see |23| for an analogy. We believe that the mechanisms 
of gauge-symmetry breakdown presented here have a lot 
to say about electroweak symmetry breaking without the 
conventional prediction of Higgs-particle excitations. 
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